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We explicitly calculate nonthermal gravitino production during the preheating period in the infla- 
tionary Universe. Contrary to earlier investigations, we consider a two-fleld model to separate the 
mechanisms of supersymmetry breaking and inflation. We show that the superpartner of the infla- 
ton is signiflcantly generated, while the gravitino production is considerably smaller. Nonthermal 
production of gravitinos seems thus less worrisome than recently claimed. 
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I. INTRODUCTION 

Inflationary models, when embedded in particle- 
physics motivated schemes, are usually constrained by 
the requirements of consistency with the phenomenol- 
ogy of the later evolution of the Universe. In particular, 
in the context of supergravity theories the parameters 
of the inflationary sector have to be restricted in order 
to avoid the thermal overproduction of gravitinos at re- 
heating. These constraints, in the case of models with 
gravity-mediated supersymmetry breaking, can be sum- 
marized as an upper limit on the reheating temperature 
of the order of 10^ GeV (for a review, see Q). In the last 
years it was shown that, during the stage of coherent os- 
cillations of the inflaton condensate right after inflation, 
preheating can lead to an efficient nonthermal production 
of fermions Q. In particular, it was recently argued |^|^] 
that the nonthermal production of 1/2-helicity compo- 
nent of gravitinos can in some cases be much more effi- 
cient than the thermal one, thus worsening the gravitino 
problem. The investigation that led to the present letter 
was actually triggered by this last statement. Explicit 
calculations of the amount of gravitinos produced at pre- 
heating were performed so far only in models with one 
superfield and without supersymmetry breaking in the 
vacuum. Therefore, the conclusions about the produc- 
tion of longitudinal gravitinos in these models could be 
misleading, since there is no longitudinal gravitino in the 
vacuum of the theory. Thus, one might wonder whether 
preheating could actually lead to a production of harm- 
less inflatinos rather than of dangerous gravitinos. 

In order to discriminate between inflatino and grav- 
itino production it is necessary to consider more realistic 
schemes. The simplest possibility is to consider two sep- 
arate sectors, one of which drives inflation, while the sec- 
ond is responsible for supersymmetry breaking today. In 
the present analysis we study the situation in which the 
two sectors communicate only gravitationally. The calcu- 
lation of nonthermal gravitino production in this scheme 
requires substantial work. First, one has to develop a new 



formalism, to be able to clearly define and compute the 
production in systems with several coupled fields. Sub- 
sequently, an extended numerical investigation has to be 
carried out to obtain reliable results. In fact, this pro- 
cedure led us to conclusions that do not confirm earlier 
analytical estimates 0. 

The paper is structured as follows. First, we define 
the specific model considered, analyzing the evolution of 
the scalar fields. Then, we discuss how the definition of 
the gravitino varies in time according to the evolution of 
the background. We finally compute the spectrum of the 
gravitinos produced at preheating, concluding that what 
in the previous works was believed to be gravitino over- 
production should rather be regarded as inflatino pro- 
duction. 



II. GRAVITINO PRODUCTION IN A TWO 
FIELD MODEL 

The system we are considering has the matter content 
of two superfields $ and S, with superpotential 
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and with minimal Kahler potential JC = <& + 5*^ S*. 
The field (f) (that is, the scalar component of $ ) acts as 
the inflaton. As it is known, the corrections from the 
Kahler potential to the scalar potential are very relevant 
when ((/)) > Mp. This is a common problem for super- 
symmetric theories of inflation, where i^-terms usually 
spoil the flatness of the inflaton potential during infla- 
tion (for a review, see [^). This is in particular true for 
the above superpotential (|^), so that additional contri- 
butions - possibly additional scalar fields (as for example 
in D term inflation [^) - must be relevant during infla- 
tion. However, we are interested in the dynamics of the 
system after inflation, when (0) < Mp, and supergrav- 
ity corrections are not important. We then assume that 
in the reheating stage only one of the fields which were 
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driving inflation is still relevant, and we refer to it as the 
inflaton. We also assume ^ 10^^ GeV, as required 
by the COBE normalization of the CMB fluctuations for 
the "usual" chaotic inflation. 

The superfleld S leads to the breaking of supersym- 
metry in the true vacuum owing to its "Polonyi" super- 
potential 0]. By imposing P — (2 — \/3) Mp, one can 
indeed break supersymmetry while retaining a vanish- 
ing cosmological constant in the true vacuum, where the 
fields s (the scalar component of the superfield S) has ex- 
pectation value of the order of Mp. The gravitino mass 
in the vacuum is of the order /Mp. In order to have 
a gravitino mass of about 100 GeV (that is the expected 
value for the gravitino mass in gravity-mediated super- 
symmetry breaking models), /i ~ 10^° GeV is required. 

The superpotential W describes a system where su- 
persymmetry is broken in a "hidden" sector and is trans- 
mitted gravitationally to a "visible" one, to whom the 
inflaton belongs. This separation between the two sec- 
tors is suggested by the two different mass scales which 
characterize inflation (inflaton mass ^ 10^"^ GeV) and su- 
persymmetry breakdown (gravitino mass ^ 10^"'^ GeV). 

Right after inflation the fleld is oscillating about the 
bottom of its potential with frequency proportional to 
The time dependent expectation value of 4> acts as 
an effective mass for the Polonyi scalar, which has vanish- 
ing expectation value at this stage. In this initial period 
the (time dependent) expectation value of (p is the main 
source of supersymmetry breaking. The amplitude of the 
oscillations of the fleld (j) eventually decreases, due to the 
expansion of the Universe, and for times of the order of 
TOj^^2 the Polonyi scalar starts rolling down towards its 
true minimum and then oscillates about it. ^ 

The system is thus governed by two time scales. At 
"early" times, of the order of m^^, the only relevant dy- 
namics is the one of the inflaton sector, that is also the 
main source of supersymmetry breaking. At "late" times, 
much larger than rn'^J^^ the system behaves as if it was 
in its true vacuum, and supersymmetry is broken by the 
Polonyi sector. To be more specific, we define the dimen- 
sionless parameter /i^ = /j,^/ {m^ Mp) ~ 1713/2 /m^, that 
gives the ratio of the two time scales in the system. If su- 
persymmetry is supposed to solve the hierarchy problem, 
p,^ should be of the order of 10~^^. Such a small param- 
eter implies a very large difference between the two time 
scales of the problem, which is a source of technical diffi- 
culties in the numerical computations. As a consequence, 



we could not study the evolution of the system for such 
a small value of /i^. Thus, we kept it as a free parameter 
and we studied how a variation of fi^ affects the scaling 
of the relevant quantities. 

To make the difference in the two time scales that gov- 
ern the system manifest, we show in fig. |^ the evolu- 
tion of the inflaton and the Polonyi scalar for the case 
= 10-2. 




FIG. 1. Evolution of the two scalar fields 
10" ^ 



We can indeed see the two different stages of the scalar 
evolution, the flrst of which is characterized by the in- 
flaton oscillations, while the second is governed by the 
dynamics of the Polonyi fleld. While (s) gives the su- 
persymmetry breaking in the vacuum, both scalar flclds 
contribute to break supersymmetry during their evolu- 
tion. In particular, both their kinetic and potential en- 
ergies contribute to the breaking, as emphasized in 
Following H], we deflne 
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dt 
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with rrii = exp (/C Mp ^/2) [d.,W + Mp ^diJCW]. The 
quantities fi give a "measure" of the size of the supersym- 
metry breaking provided by the F term associated with 
the z— th scalar fleld. More precisely, we will be interested 
in the normalized quantities = ff/ (/f -I- /|) , which 
indicate the relative contribution of the two scalars. 

In flg. 1^ we show the evolution of and r^. for the 
specific case fi^ = 10~ ^ . As it occurs for all values of 
A^ , in the initial stages ~ 1 , while rj ~ 1 at the end. 
The regime of equal contribution is around t = A~ ^ • 



*In general, the presence of the Polonyi sector will affect the 
inflationary dynamics (as discussed in Again, since we 

are not interested in the inflationary expansion itself, we will 
not consider this issue here. 

^ There is of course a possible moduli problem associated 
with these oscillations. However, we do not consider this issue 
here. 
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FIG. 2. Relative contribution of the two scalar fields (j) and 
s to the supersymmetry breaking during their evolution. As 
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Let us now consider the fermionic content. We denote 
the fermions of the two chiral multiplets by (the "in- 
flatino") and s (the "Polonyino"). One hnear combina- 
tion of them is the goldstino v , while the one orthogonal 
to V is denoted by T . Initially, v = (j) , while u = s at 
late times. In addition, we have the gravitino field, whose 
longitudinal and transverse component are denoted by 
and by iff , respectively. The transverse component is 
decoupled from the other fermion fields, and its quanta 
are produced only gravitationally P,pp^. It will not 
be considered in the remainder of this work. The lon- 
gitudinal gravitino component, which in the super-higgs 
mechanism is provided by the goldstino v , is however 
coupled with T . Equations of motion for the coupled 
system are given in H. The step from these equations 
to the occupation numbers of the fermions 9 and T is far 
from trivial and requires a significant extension of the 
existing formalism for nonperturbative production in the 
one field case. We will present the details in a separate 
publication [p^ . In the remainder of this work we will 
only present our results. 
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FIG. 3. Evolution of the masses of the two fermionic eigen- 
states. As in fig. |l|, /l^ = 10~ ^ . Notice the different normal- 
izations for the two masses. 



As a starting point, one has to diagonalize (at each 
time) the coupled 9-T system. We denote the two 
fermionic mass eigenstates by ipi ^-nd tl'2 ■ In fig- H we 
show [] the evolution of their masses for the specific case 
fi? ~ 10^ ^ . The most relevant information which can be 
derived by this evolution is very clear: at late times the 
fields ipi and il)2 have, respectively, the mass of the in- 
flatino and of the gravitino field. That is, at late times we 
have the identification ■0i=0 = T,^2 = ^ {s — v — Q , 
being the goldstino). This situation is orthogonal to the 
initial one, when ipi = 9 , il)2 = s = T {(f) — v — 0). 
At intermediate times, when supersymmetry is broken by 
both scalar fields, the gravitino is a mixture of ipi and ip2 ■ 
To qualitatively appreciate the evolution of the gravitino 
occupation number, we may consider Ng = ri Ni +r2 N2 
and the orthogonal combination A^x = r2 Ni + ri N2 , 
where are the relative contributions of the two scalars 
to supersymmetry breaking defined above. 




FIG. 4. Evolution of Ng and A^t for /i^ = 10"^ and 
k = . See the text for details. 

We show in fig. ^ the evolution of Ni , N2 , Ng , Nr 
for modes of comoving momentum k = and for fl^ = 
10"^ . Notice that (by construction) Ng = Ni at early 
times, while Ng = N2 at late ones. We also see that ipi is 
populated on time scales ^ , while '02 on time scales 

A~ ^ ^ - This feature is common for all fi^ . We 
remark that the identification 9 = ri'ijji+r2 ^"2 should be 
taken only as a qualitative indication. However, the most 
relevant identification = ■02 at late times is a rigorous 
one, as should be clear from the above discussion. 

We are now ready to present our most important re- 
sult: the occupation number of T and 9 at the end of the 
process. We show them in figs. || and ^, respectively. || 



"^From fig. ^ one may be tempted to identify tpi = (f> and 
i/>2 = S . Although this identification is rigorous only at the 
beginning and at the end of the evolution, it can be used for 
an "intuitive" understanding of the system. 



^These spectra are shown at the time t — 10 p, m. 
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The time required for the numerical computation in- 
creases hnearly with /i^ , so that we could perform it up 
to jl = 10"^ . The realistic case fi^ — 10"^^ is far from 
our available resources, so we kept /i^ as a free parameter. 
The case fi^ — 10~ can be clearly extrapolated from 
the results we are going to discuss. Moreover, the case 
/t^ = can be studied analytically, and it agrees with the 
limit jl ^ Q that one deduces from the numerical results. 
For /i = , only the inflatino is produced. The mass of 
the Polonyi fermion does instead vanish identically, so no 
quanta of this particle are produced at preheating [ p^ . 
Notice that [i — Q corresponds to a situation with un- 
broken susy in the vacuum, and it reproduces the models 
with one single field studied so far. We stress that in this 
case only the inflatino is produced at preheating. 
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FIG. 5. Spectrum of inflatinos at late times. 
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FIG. 6. Spectrum of gravitinos at late times. 



the /i = 10"^ - 10"" cases we have continued the evolution 
further, until the spectra stop evolving. We have found that 
the spectra shown in fig. |^ coincide with the final ones, while 
Ne very slightly decreases for t > 10/1" ^ ■ Thus, the 
results shown in fig. ^ give an accurate upper bound on the 
final gravitino abundance. 



In fig. K| the spectrum of ^"1 in shown for jj? = 
10"^ 10"^ 10"'', 10-^ 10"^. It is apparent that the 
main features of the spectrum are independent of the 
value of . The reason for this is that this eigenstate 
is associated to the inflatino, that is produced by the co- 
herent oscillations of the inflaton. The inflaton dynamics 
occurs on time scales of the order of m^^, and is indepen- 
dent on ji^ . The spectrum given in fig. ^ shows a more 
significant dependence on ji^ . In particular, as in the 
case of the spectrum in fig. |^, the occupation number is 
of the order of unity if the comoving momenta are smaller 
than some cut-off fc* . At variance with the spectrum of 
V'l , however, here the cut-off depends on the value of jx^ . 
The total number of particles produced in this case is a 
increasing function of /i^. 

It is possible to estimate in a very simple way how fc* 
scales with fi^ . The cut-off in the spectrum (in terms of 
the comoving momentum) of particles produced at pre- 
heating is indeed generically proportional to ^/iha, where 
m is the time-derivative of the effective mass of the pro- 
duced particle, and a is the scale factor of the Universe at 
the time of preheating. In the case of the Polonyino, the 
typical mass scale is rn3/2, while the typical time scale is 
mgyj- As a consequence, ^frh ~ m3/2- When the energy 
in the Universe is dominated by the sinusoidal oscillations 

of a massive scalar field, a (t) ^ t^/^ ^ {^11^/2) Col- 

1/3 

lecting all these estimates, we obtain ~ (/i^) 
Although this is quite a simple estimate, it shows very 
good agreement with the numerical results. Moreover, it 
agrees with the fact that iV2 = for p,^ — 0. 

In conclusion, our calculation confirms that one 
fermionic eigenstate is efficiently produced at preheating. 
However, as we have initially remarked, in more realistic 
models (i.e. with more than one chiral superfield and su- 
persymmetry broken in the vacuum) this fermionic field 
should be regarded as the "inflatino" rather than the 
"gravitino" . We have shown that gravitino production is 
significantly reduced if the sector responsible for super- 
symmetry breaking today is coupled only gravitationally 
to the one responsible for inflation. The reason for this 
small production is related to the smallness of the super- 
symmetry breaking scale (i.e. of fi^) with respect to the 
scale of inflation. 
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